Pagliarani et al. BMC Plant Biology 2013, 13:51 
http://www.biomedcentral.com/1471-2229/13/51 



Plant Biology 



METHODOLOGY ARTICLE Open Access 



A qRT-PCR assay for the expression of all Mai d 7 
isoallergen genes 

Giulia Pagliarani 1,2 * + , Roberta Paris 2,41 , Paul Arens 1 , Stefano Tartarini 2 , Giampaolo Ricci 3 , Marinus MJ Smulders 1 
and W Eric van de Weg 1 



Abstract 

Background: A considerable number of individuals suffer from oral allergy syndrome (OAS) to apple, resulting in 
the avoidance of apple consumption. Apple cultivars differ greatly in their allergenic properties, but knowledge of 
the causes for such differences is incomplete. Mai d 1 is considered the major apple allergen. For Mai d 1, a wide 
range of isoallergens and variants exist, and they are encoded by a large gene family. To identify the specific 
proteins/genes that are potentially involved in the allergy, we developed a PCR assay to monitor the expression of 
each individual Mai d 1 gene. Gene-specific primer pairs were designed for the exploitation of sequence differences 
among Mai d 1 genes. The specificity of these primers was validated using both in silico and in vitro techniques. 
Subsequently, this assay was applied to the peel and flesh of fruits from the two cultivars 'Fiorina' and 'Gala'. 

Results: We successfully developed gene-specific primer pairs for each of the 31 Mai d 1 genes and incorporated 
them into a qRT-PCR assay. The results from the application of the assay showed that 1 1 genes were not expressed 
in fruit. In addition, differential expression was observed among the Mai d 1 genes that were expressed in the fruit. 
Moreover, the expression levels were tissue and cultivar dependent. 

Conclusion: The assay developed in this study facilitated the first characterisation of the expression levels of all 
known Mai d 1 genes in a gene-specific manner. Using this assay on different fruit tissues and cultivars, we 
obtained knowledge concerning gene relevance in allergenicity. This study provides new perspectives for research 
on both plant breeding and immunotherapy. 
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Background Of the 4 classes of allergens currently identified in 
Apple allergy is an issue for a growing number of apple, Mai d 1 is thought to be the major allergen in 
European citizens. As one of the most prevalent food Central and Northern Europe [4,5]. At the genetic level, 
allergies, apples rank fourth out of 24 foods examined Mai d 1 is a complex gene family composed of 31 differ- 
in an extensive Pan-European survey and first among ent loci, each of which codes for a different isoallergen 
Rosaceae fruits [1]. Thus, although the apple is generally [6]. Moreover, for each isoallergen gene, there are a 
a healthy component in the human diet, an increasing series of slightly different alleles that might encode for 
number of individuals cannot eat this fruit. The allergenic isoallergen variants, which increases the variability in 
properties of apple cultivars differ greatly [2,3], but Mai d 1 proteins [6-8]. Accumulating evidence has 
knowledge of the genetic basis for low and high allergen- shown that isoallergens might differ greatly in their 
icity remains incomplete. allergenic properties, but it is still unclear which of these 

proteins are more involved in allergy. Several approaches 

have been used to quantify Mai d 1 content or gene 
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high temperature and proteases hinders its proteomic 
analysis. Moreover, the food matrix and contaminants 
might affect the protein extraction. Until now, proteomics 
have been primarily used to quantify the total amounts of 
Mai d 1 content in apple fruit, without distinguishing 
isoallergens or variants or making distinctions within an 
incomplete pool of isoallergens [10-14]. Currently, 
there are only a few recombinant allergens derived 
from fruits and vegetables that are commercially available 
for immunological detection [15], with variable antibody 
specificities [16]. 

PCR-based expression studies are not subject to these 
limitations, and in particular, quantitative Real-Time PCR 
(qRT-PCR) is a fast, highly sensitive and reproducible tech- 
nique to study gene expression. Previous studies of some 
Mai d 1 genes revealed the tissue- and cultivar-specific 
expression of Mai d 1 genes [17-19], and differential effects 
of environmental conditions [20,21] on the transcription of 
these genes have been reported. However, these studies 
have not covered the entire gene family, nor have they 
sufficiently demonstrated the gene specificity of the 
PCR primer pairs used. 

Thus, there is a need to characterise the role of each 
individual isoallergen to understand the apple allergy 
mechanism. Therefore, the aim of this study was to 
examine the expression of these genes by implementing 
the qRT-PCR approach. A comprehensive, robust, sensitive 
and affordable assay for studying the expression of all 
31 known Mai d 1 genes individually was developed. 
We successfully used this assay to generate the 
complete expression profile of all Mai d 1 isoallergen 
genes in the fruits of two cultivars. 

Results 

Alignment of Mai d 1 genes 

A total of 380 Mai d 1 DNA and EST sequences were re- 
trieved from the literature and databases (Additional file 1) 
covering all 31 known Mai d 1 genes described in [6]. 
Many sequences obtained from different apple cultivars 
were already annotated as Mai d 1 alleles. For the others, 
according to the level of similarity, it was possible to 
identify new alleles of known genes. Subsequent to align- 
ment, the level of similarity among the coding sequences 
(cds) of different genes ranged from 53.1% to 97.7% and 
from 95% to 99.8% for the different alleles within a gene. 
At the protein level, the sequence identity between 
isoallergens ranged from 37% (Mai d 1.08 and Mai d 
1.1 1A, 102 AA substitutions) to 96% (Mai d 1.06A and 
Mai d 1.06D, 7 AA substitutions). The alignment of the 31 
coding sequences retrieved from the 'Golden Delicious' 
whole genome sequence (Additional file 2) was used to 
generate a phenetic tree (Figure 1). This tree showed 5 
clades, 4 of which have been previously described and 
characterised as subfamilies I-IV [7] and one clade, 
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Figure 1 Mai d 7 phenetic tree. Neighbour-joining tree based on 
the coding sequence of the 31 Mai d 1 isoallergen genes retrieved 
from the 'Golden Delicious' genome sequence (Additional file 2). 
The sequences are presented using the names of the related genes 
and the accession numbers obtained from Apple GBrowse - Malus x 
domestica vl.O [22]. The Roman numerals (l-V) and colored lines 
identify the subfamilies. 



subfamily V, that has been classified for the first time 
in the present study. The three genes in subfamily V 
(Mai d 1.11A, Mai d 1.11B and Mai d 1.12) are the 
most distant within the Mai d 1 family. 

Primer design 

The development of this assay began with the design of 
specific primer pairs. The alignment of the 31 Mai d 1 
sequences from 'Golden Delicious' (Figure 2A) highlights 
SNPs specific for only one gene. The robustness of these 
SNPs across alleles of the same gene was examined 
among the other allelic sequences in the alignment, 
using the 'Golden Delicious' sequences as a framework. 
An average of four gene-differentiating SNPs per sequence 
was detected and only these SNPs were exploited for 
primer design. For example, the reverse primer for Mai d 
1.02 is located in a region in which a SNP differentiates the 
Mai d 1.02 gene from all other genes (Figure 2A). The 
alignment of all allelic sequences for Mai d 1.02 (Figure 2B) 
shows that no additional allele-differentiating SNPs were 
present in this region; thus, this region was an excellent 
candidate for primer development. Figure 2B shows an 
example of the reverse primer for Mai d 1.01. In addition 
to the gene-differentiating SNP (SNP nr. 2) targeted 
for primer design, this primer contained one allele- 
differentiating SNP (SNP nr. 6) at the 5' end. In general, 
regions containing only gene-differentiating SNPs were 
preferred for primer design; however, these regions were 
not always available. Thus, 12 of the 31 primer pairs 
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Figure 2 Primers design strategy. A) Part of the coding sequence alignment of the 31 Mai d 7 isoallergen genes was retrieved from the 
'Golden Delicious' genome sequence (Additional file 2). The sequences were named according to their related Mai d 1 genes and the accession 
numbers retrieved from the Apple GBrowse - Malus x domestica vl.O [22]. The white marked nucleotides highlight the mismatched residues in 
the consensus sequence. The SNPs specific for the loci Mai d 1.01 and -7.02 are indicated with red circles. B) Part of the alignment on all allelic 
gDNA sequences for Mai d 1.01 and -7.02 was obtained from the literature [8]. The sequences of the two reverse primers for Mai d 1 .01 and -7.02 
are also included in the alignment. The sequences were named according to their related Mai d 1 genes, ID numbers from the database and 
genotypic origins: FS: Fiesta; RD: Red Delicious; DS: Discovery; PM: Prima; GD: Golden Delicious; JO: Jonagold; FJ: Fuji; IM: Ingrid Marie; SC: 
Suncrisp; RG: Royal Gala. The SNP positions are indicated with successive numbers from 1 to 1 1. The red boxes indicate the two locus- 
differentiating SNPs exploited for primer design; the green box highlights an allele-differentiating SNP. 



contained allele-differentiating SNPs (Table 1). To ensure 
that these SNPs would not affect the PCR amplification, 
the allele-differentiating SNPs were accepted only if 
positioned at the 5' end of only one of the two primers. 
The procedure for the primer selection was performed 
for each of the 31 Mai d 1 genes (Table 1). 

Primer validation and qRT-PCR optimisation 

The gene-specificity of the designed primer pairs was 
validated in four ways. Firstly, an in silico validation 
through blasting the primer sequences to known reference 
sequences ensured that these primers perfectly matched 
only with sequences corresponding to the target genes 
(data not shown). Only the primers that generated perfect 



matches were assessed in the second validation step 
performed through end-point PCR on genomic DNA; 
only the primer pairs that produced a single clear band 
were maintained. Thirdly, the direct sequencing of the 
Mai d 1 amplicons obtained from 10 apple cultivars was 
performed to ensure that only the specific target 
sequences were amplified (Additional file 3). Amplicons 
generated from gDNA were sequenced because of the 
higher level of complexity of gDNA compared to cDNA 
and this allowed to guarantee the primers specificity for 
all the Mai d 1 members, independently from their 
expression. Only the primers that had a unique sequence 
corresponding to the target gene were assessed in the last 
step of validation. In some cases, this specificity was 
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Table 1 Mai d 7 specific primer pairs for qRT-PCR analysis 



Gene 


Primer name 


Sequence 5 -3 


SNP 
cons 3 


SNP 
gene b 


SNP 
allele" 


Start 
position 


Length 
(bp) 


Primer 
cone. 


Ta 


Tm 


Slope 


Maid 1.01 


qMd1.01/02F 
qMd1.01R 


GATOAAGGAGATGCTITGACA 
GTAATGACTGATGCTCTOATGG 


5 


1 


1 


258 


103 


100 


63 


80.5 


-0.079 


Mai d 1.02 


qMd1.01/02F 
qMd1.02R 


GATOAAGGAGATGCTITGACA 
TOGTGTGGTAGTGGCTGAT>A 


4 

1 


1 


- 


258 


111 


100 


62 


81.0 


-0.206 


Mai d 1.03A 


qMd1.03AF 
qMd1.03AR 


ATCrGAGTOACCTCCGTCATr 
ACTGCTTGTGGTGGAATCTTT 




2 
1 


1 


21 


96 


70 


63 


81.0 


-0.057 


Mai d 1.03B 


qMd1.03BF 
qMd1.03BR 


tgi i i i cacatacgaatccga/4 
tgatcto™atggtoctacgc 


1 


1 
1 


1 


6 


167 


100 


63 


83.5 


-0.570 


Mai d 1.03C 


qMd1.03CF 
qMd1.03CR 


CTCCGAAACAATOAGAAAATCTG 
GCTGGTGCTCTOATGATGC 


1 


1 


1 
- 


276 


79 


100 


63 


80.5 


-0.120 


Mai d 1.03D 


qMdl .03D/EF 
qMd1.03DR 


ATACGAATCCGAGTCACCTCT 
ATCITCrTAATGG^CCAACTCCr 




1 




15 


156 


70 


62 


83.0 


0.005 


Mai d 1.03E 


qMdl .03D/EF 
qMdl .03ER 


atacgaatccgagtoacctct 
tcaccgaagtoatctc^aa™ 




1 


2 
- 


15 


169 


70 


62 


83.0 


-0.116 


Mai d 1.03F 


qMd1.03FF 
qMd1.03FR 


CACAGAATOACGGGGTG 
CCGGAAGCGACCAACTTA 


2 
2 


2 
1 


1 


208 


119 


70 


63 


81.0 


-0.122 


Mai d 1.03G 


qMd1.03GF 
qMd1.03GR 


ATTATCAAGAGCACCAGrCACTACr 
TCCAAGAGGTAGTOTCAATCAA 


2 
1 


2 


- 


337 


122 


70 


62 


81.0 


-0.254 


Mai d 1.03H 


qMdl .03HF 
qMd1.03HR 


AAAATOX5CTACGAGACTAAGTTGA 
TGGTGCTCCAAGAGGTAGTT7 


3 
1 


2 

1 


- 


277 


173 


100 


61 


83.5 


-0.431 


Mai d 1.031 


qMd1.03IF 
qMd1.03IR 


CCCCAAGATOCACCACA7 
GCCACCAACTOGTCTCGTAACAyA 


2 


1 


- 


93 


228 


100 


61 


81.5 


-0.230 


Mai d 1.03J 


qMd1.03JF 
qMd1.03JR 


GCATCACCCACTACCACACvA 
CGAGCTGTAGG4GTC7TGGT7 


2 
3 




- 


347 


134 


70 


61 


82.0 


-0.105 


Mai d 1.03K 


qMdl .03KF 
qMdl .03KR 


CATCAGCCACTACCACAGAAA 
TGTATGCATCCTGGTGCTC7 


2 
2 


1 


1 


348 


128 


100 


61 


81.5 


-0.431 


Mai d 1.04 


qMd1.04F 
qMd1.04R 


GGGTATGTTAAGCMAGGGTOl 

TGATCTCAACATCACCCTTAGC 


5 
1 


1 




196 


193 


100 


61 


80.5 


-0.103 


Mai d 1.05 


qMd1.05F 
qMd1.05R 


ATCAAACCACTAGTCACTGCCAT 
GGTTGGCCACAAGGTAGGTr 


4 
6 




1 
- 


343 


124 


70 


63 


82.5 


-0.141 


Mai d 1.06A 


qMdl .06AF 
qMdl .06AR 


CTATAGCTATAGCTOATOAAGGG 
rrCCAACCTOACATGTTCTOT 


5 
3 




1 


243 


167 


100 


61 


80.5 


-0.203 


Mai d 1.06B 


qMd1.06BF 
qMd1.06BR 


AAACCGAATACGCATCCdTr 
ACAGTmGACTGCTOTGGAG 


3 
6 


3 
- 


1 
- 


20 


106 


100 


61 


81.5 


-0.012 


Mai d 1.06C 


qMald1.06CF 
qMald1.06CR 


GCTCCACAAGCAGTCAAAACT 
TCAACCTTGTGCTOACATA>ACT^ 


5 
3 


2 


- 


103 


116 


70 


63 


80.5 


-0.250 


Mai d 1.06D 


qMd1.06DF 
qMd1.06DR 


CCCTCCTGCTAGGTOTA7T 
TCCCTCGAGAAmCAACAG 


2 
6 


2 




42 


100 


70 


61 


80.5 


-0.005 


Mai d 1.07 


qMd1.07F 
qMd1.07R 


CAACmGTGTACCAGTACAGTGTC 
TAGTGGCTGATGCTC1TGAT>AAC 


2 
2 


2 
1 




234 


126 


100 


61 


81.5 


-0.201 


Mai d 1.08 


qMd1.08F 
qMd1.08R 


7CTOGGTGAAGGTAGCACA>A 
ACCC™GTGTGGTAGTGGC^\r 


1 


2 
2 




173 


200 


100 


61 


81.5 


-0.390 


Mai d 1.09 


qMd1.09F 
qMd1.09R 


^CACATACGAATCCGAGTC 
GGATCTCAACGCTCTCAC>4 


2 


1 
1 




8 


126 


100 


61 


84.0 


-0.265 
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Table 1 Mai d 7 specific primer pairs for qRT-PCR analysis (Continued) 



61 83.0 -0.158 



62 79.5 



61 78.5 



-0.018 



Maid 1.10. qMaldl.lOF CAAGGCmCATCCACGAC 5 2 1 60 157 100 

qMaldl.lOR GATTCTGTGOTrACAyAACCCT 4 3 

Maid 1.11 A qMald 1 .11 AF GGAGGATGCATCTGTCATTTG 11 1 - 343 1 30 1 00 

qMaldl.HAR CCATGAG>ATAGGC1TCCAA>AACT 8 2 

Maldl.llB qMdl.HBF CAGCACATACAAAGCCAAAGAC 8 1 - 363 125 100 

qMdl.HBR TITATGCGCGAGGGTGrG 6 1 

Maid 1.12 qMd1.12F GCmCACnTTGGlTGAAGGAGMC 4 3 247 171 100 

qMd1.12R CaGCCAGCmrATTTCTTCC 5 4 

Maid 1.1 3A qMd1.13AF GTGTOGAACCATCAAGAAGAT7AG 2 2 - 149 124 100 

qMd1.13AR ACATCTCC1TCAATCAAACTGTMT 1 1 

Mald1.13B qMd1.13BF CGAAGATAACmGTCTACAACCAT 2 1 - 258 137 70 

qMd1.13BR GCTCTOOTGATCTCAACATCT7 1 1 

Maid 1.1 3C qMd1.13CF GAATTCGCCTCAGTCTCC^l 5 1 - 25 186 70 

qMd1.13CR GTGCTTCACATAGCTGTATTCAGT 3 1 

Maid 1.1 3D qMd1.13DF TGTOGAACCATCAAGAAGATAAGT 2 1 - 150 124 100 

qMd1.13DR GACATCTCGTCAATCAAACTGTAG 2 1 

Maid 1.14 qMdl.HF GGTGAAGGGAGTGAATACAACTATA 6 1 - 178 185 100 

qMdl.HR TGGTAATGGCTAATGTTOTGATAC 2 2 2 

a SNPs to the consensus sequence showed in bold; b gene-differentiating SNPs showed in italics and c allele-differentiating SNPs showed by underlining. The last 

column contains the slopes of the curves obtained by plotting log input vs ACt (Ct Mai d 1 gene - Ct actin). The values between -0.1 and 0.1 are underlined. 



61 81.0 -0.106 



62 76.0 -0.227 



61 78.0 -0.216 



61 81.5 0.003 



61 82.0 0.087 



-0.163 



61 79.0 -0.470 



obtained after increasing the annealing temperature 
(Ta) within the range of 61-63°C and/or decreasing the 
final primer concentration from 100 nM to 70 nM. 
Some of the amplicon sequences revealed new Mai d 1 
alleles for the known genes (Additional file 3). Fourthly, 
the amplification specificity was validated through 
melting curves obtained with a single high peak (Figure 3), 
indicating the absence of nonspecific amplifications 
or primer dimers. The primer pairs that met this last 
criterion were selected for the qRT-PCR assay. The 
melting temperatures (Tm) for these primers are reported 
in Table 1. 

Expectedly, significant variations in PCR efficiency 
were detected among the Mai d 1 amplicons (Table 1, 
slopes in the last column) due to the forced positions of the 
primers to the targeted SNPs and variations in the PCR 
conditions to ensure specificity. However, the "Standard 
curve method" was chosen to neutralise these differences in 
efficiency. Particularly, this method has to be used when 
the slope values of the curves obtained by plotting log input 
vs dCt (Ct Mai d 1 gene - Ct actin) are outside the -0.1 to 
0.1 range [23], as in the case of this study (Table 1). The 
standard curves were of high quality, as a linear relationship 
between the input DNA and the Ct values across the 
standard samples (serial dilutions) was observed, with a 
squared regression coefficient close to 1 for all genes. 
Moreover, the use of standard samples in all the qRT-PCR 
reactions facilitated the evaluation of the reproducibility 
among experiments and the integration of the data. 



Finally, the comparison of the gene expression in 
'Fiorina and 'Gala of actin, UBC and GAPDH was 
performed in order to choose the most suitable reference 
gene for this study. Actin showed a highly stable expression 
among tissues and cultivars (Figure 4) and thus, it was an 
optimal reference gene for this study. 

Application: Mai d 7 expression profiles in different apple 
tissues 

In this study, we applied this new qRT-PCR assay to the 
peel and flesh of apple fruits from the cultivars 'Fiorina 
and 'Gala. These two cultivars were selected because 
they differ in their allergenic characteristics in skin prick 
tests [24]. Moreover, both cultivars have been extensively 
used in the breeding of new apple cultivars and are 
subject to GM studies concerning the enhancement of 
'Gala resistance to the main disease in apple production, 
the fungus Venturia inaequalis, through gene transfer 
without raising Mai d 1 levels [25]. 

Among the 31 Mai d 1 genes, 11 were not expressed 
in any fruit tissue. These genes included all of the genes 
of subfamily II and several genes of subfamilies III and 
IV (respectively Mai d 1.06D, -1.13C and -1.14, and Mai 
d 1.03B, -1.03H, -1.031, -1.03], and -1.09) (Figure 5). The 
20 genes expressed in these fruits included all of the 
genes of subfamilies I and V and some of the genes of sub- 
families III and IV (respectively Mai d 1.06A, -1.06B, - 
1.06C, -1.13A, -1.13B, -1.13D and -1.03A, -1.03C, -1.03D, - 
1.03E, -1.03K, -1.03F, -1.03G, -1.07, -1.08). A significant 
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variation in the transcript level among genes, tissues and 
cultivars was observed (P<0.0001). 

Concerning the variations in expression among the 
genes, the two genes of subfamily I (Mai d 1.01 and -1.02) 
showed the highest expression level, which were 10 to 
10,000 times higher than any other Mai d 1 gene, and the 
combined total amount of RNA for these genes exceeded 
the amount of all other genes combined. Mai d 1.01 was 
the most expressed gene, with a normalised expression of 
approximately 60 A.U. in the fruit peel of 'Fiorina . Mai d 
1.03C was least expressed gene, with a transcript level 
close to the limit of detection. 

The variations in gene expression between fruit tissues 
were evident, with the expression in the apple peel being 
significantly higher than in the flesh for both cultivars. 
The genes of subfamily V were the only exception, as 



Mai d 1.11 A and -1.1 IB showed similar levels of expression 
in the peel and flesh, and Mai d 1.12 was more expressed 
in the flesh than in the peel in 'Gala fruits. 

A comparison of the two cultivars revealed that 
'Fiorina generally shows a higher level of expression 
than 'Gala', which was consistent with a previous study 
on the expression of a limited number of Mai d genes [19]. 
Specifically, for 15 out of 20 genes, the expression was 
higher in 'Fiorina! with an average increase of 60%. The 
remaining 5 genes, Mai d 1.03A, -1.03K, -1.06A, -1.06C 
and -1.07, belonging to subfamilies III and IV, were more 
expressed in 'Gala than in 'Fiorina . 

Discussion 

The results obtained from the assay described in this paper 
provided insight into the role of individual isoallergens in 
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Figure 4 Ct values of three putative reference genes in 'Fiorina' 
and 'Gala' peel and flesh. The Ct values obtained after qRT-PCR 
amplification with the primers for actin, UBC and GAPDH in 'Fiorina' 
and 'Gala' samples are plotted in the chart. Each data point is the 
average of two biological replicates and the bars indicate the 
standard deviation. 



Mai d 1 -related apple allergy. The assay was based on a 
complete set of primers that were suitable for expression 
studies for each Mai d 1 isoallergen gene using qRT-PCR. 
All primers passed several assessments of validation and 
optimisation and can now be used in a broad range of ex- 
periments. This approach can be relatively quickly adopted 
to other crops whose allergen genes have been mapped, 
a reference genome sequence is available and some 
knowledge of the genes and allelic diversities is available, 
such as for the peach Pru p 1 genes [6,22,26]. 

The qRT-PCR assay 

The study of gene expression is an informative approach, 
and the first application of this assay confirmed Mai d 1 
as a heterogeneous family of genes that show distinct 
expression patterns despite having highly similar sequences. 
For single copy genes, the qRT-PCR technique is a highly 
sensitive and reproducible technique, with a large dynamic 
range compared with microarray approaches [23], without 
requiring a strong knowledge of bioinformatics, expensive 
equipment or particular expertise. qRT-PCR is also 
relatively cheap, particularly when SYBR green chemistry is 
used, as in this study, initial accurate primers are selected 
and an appropriate methodology is employed. Here we 
demonstrated the suitability of qRT-PCR for the assessment 
of a large gene family, whereby we obtained gene specificity 
through the use of sequence polymorphisms at the 
3'end of at least one of the two primers and stringent 
PCR conditions. 

The calibration method can seriously affect the results 
of qRT-PCR based assays. The "Standard curve method" 
accounts for variations in PCR efficiency among primer 
pairs [23], thus facilitating the comparison and integration 
of data from different primer pairs and experiments. The 



only disadvantage is the need to include a series of standard 
samples in each qRT-PCR experiment for all the tested 
genes, including reference gene, thereby slightly increasing 
cost and time compared with other methods. 

Clarifying the role of individual genes 

Small changes in protein sequences might result in large 
differences in allergenic characteristics. A single amino acid 
change at a critical position in the epitope formation might 
completely alter the allergenicity of Mai d 1 proteins 
[27-29]. For Bet v 1, a difference of 7 amino acids between 
two isoallergens resulted in large differences in their 
allergenic and immunogenic properties: Bet v 1.0101 acted 
as a strong sensitiser, whereas Bet v 1.0401 was clearly 
associated with weaker IgE responses [30]. A similar 
situation likely exists in the apple, implying the existence 
of naturally hypoallergenic Mai d 1 isoallergens that 
induce no or mild IgE responses in allergic individuals. 
Currently, the allergenicity of single Mai d 1 isoallergens 
and their variants remains undetermined. Previous studies 
using gene expression approaches to obtain information on 
the differential expression within this gene family have 
characterised Mai d 1.02 as the most expressed isoallergen 
gene [17-20]. However, our assay is a crucial step towards 
the identification of specific candidate genes for allergen- 
icity within the Mai d 1 family through an assessment of 
the presence and unequal dispersion of all the Mai d 1 gene 
transcripts in apple fruits, thereby correlating differences in 
allergenicity among cultivars with differences in the 
expression of individual genes. 

Indeed, the results obtained from the first application 
of this assay showed that 11 of the 31 Mai d 1 genes were 
not expressed in 'Fiorina and 'Gala fruit. Therefore, the 
genes encoding these isoallergens are not involved in the 
allergenicity of 'Fiorina and 'Gala fruits and reasonably 
in the allergenicity of other cultivars, reducing the 
number of candidates for further assessments. This 
result is consistent with studies conducted in birch, 
where the Bet v 1 genes from only two of 5 subfamilies 
were expressed in pollen [31]. 

Among the 20 isoallergen genes that were expressed, 
we observed a large variation in expression between 
the tissues, genes and cultivars. We detected a higher 
expression of Mai d 1 genes in the peel than in the flesh. 
Therefore, it is reasonable to assume that peeling of apple 
fruits would remove most of the Mai d 1 proteins. Never- 
theless, peeling was helpful to a small portion of the apple 
allergic population [32], indicating that we cannot a priori 
exclude isoallergen genes based solely on their levels of ex- 
pression. Indeed, a small amount of a highly allergenic iso- 
form might provoke allergy due to high immune reactivity, 
whereas large amounts of low-immune-reactive or non- 
immune-reactive isoallergens might provoke no allergic 
responses [14]. Thus, isoallergens, such as Mai d 1.1 IB, 
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(See figure on previous page.) 

Figure 5 Expression profiles of the Mai d 7 genes in the peel and flesh of apple fruits. The genes were grouped into subfamilies and 
reported according to the order obtained in the phenetic tree. The branches of the phenetic trees for each subfamily (l-V) are reported on the 
left. The underlined genes were not expressed (charts not reported: subfamily II was not reported because none of its genes were expressed). 
The green bars indicate the results for 'Fiorina', and the orange bars indicate the results for 'Gala'. The expression levels were normalised in 
respect to octin and reported in Arbitrary Units (A.U.). Each data point is the average of the final expression values of two biological replicates. 



which are equally expressed in the peel and flesh, should 
receive full attention, despite being expressed approxi- 
mately 300 times less than Mai d 1.01. However, for in- 
dividuals for which peeling fruit is helpful, we can 
further limit the number of potentially involved genes 
to those that are only expressed in the peel 

Apple cultivars differ in allergenicity according to oral 
provocations and skin prick tests [2,3]. The current 
assay improved our understanding of these differences, 
facilitating the correlation of expression profiles across 
Mai d 1 genes to allergy responses across cultivars. 
'Fiorina exhibited lower skin prick test responses than 
'Gala in prick-to-prick tests on whole fruits [24]. Most 
of the Mai d 1 isoallergen genes were more expressed 
in the less allergenic cultivar 'Fiorina', suggesting that 
these genes play a minor role in allergy, whereas genes 
that are more expressed in 'Gala could receive priority for 
further analysis (i.e., Mai d 1.06A, - 1.06C, -1.03A, -1.03K 
and - 1.07). The indication for Mai 1.06A is of special 
interest, as a possible role for this isoallergen has been 
previously suggested [8] based on a correlation between 
the allelic composition of Mai d 1.06A and the level of 
allergenicity of a small set of apple cultivars. Extending 
the application of our qRT-PCR assay to a broader set of 
cultivars with known allergenicities will provide further 
insight into this phenomenon. We are currently applying 
this approach to apple fruits in which the expression of a 
range of Mai d 1 genes is supposedly silent [33,34] and 
for which the data from oral provocations are available 
(unpublished observations). Even if there is a general good 
correlation between proteins and transcripts abundance 
[35], quantitative information about protein content would 
be highly desirable. It is always difficult to analyse specific 
protein isoforms in the case of protein families, in particu- 
lar for Mai d 1 because of their sensitive 3D structure [9]. 
Despite that, thanks to this work it will be possible to focus 
the protein studies only on the interesting proteins, trying 
to develop Mai d 1 isoallergens-specific IgE. 

Immunotherapy can be effective for respiratory-based 
allergies, but studies have presented contradictory results 
for food allergies. The apple allergy cross-reactivity 
between pollen and foods has been exploited, but the 
effectiveness of pollen allergen immunotherapy on allergies 
to cross-reactive foods has not been confirmed [36]. Oral 
desensitisation treatment has also been examined for the 
first time, with promising results, using a mixture of 
fresh fruits from different cultivars [37]. Our assay might 



facilitate the further development of immunotherapies, as 
the allergy causing isoallergens could be identified at the 
individual level, thus leading to personalised proteins 
for use in immunotherapy. Beyond immunotherapy, 
our assay might also provide information for personalised 
diagnostics and recommendations on the safe consumption 
of specific cultivars. 

Regulation of gene expression and the biological roles of 
Mai d 1 

The Mai d 1 genes are grouped into subfamilies based on 
similarities in their genomic and amino acid sequences [6]. 
This grouping is consistent with similarities in the expres- 
sion profiles of these genes. Similar expression profiles were 
observed: i) all the genes of subfamily II were not expressed 
in fruit; ii) the two genes of subfamily I were both highly 
expressed and iii) the genes of subfamily V were equally 
expressed in the peel and flesh, suggesting sequence 
similarities in the promoter regions of these genes. 

The biological role of Mai d 1 proteins in plants remains 
unclear. These proteins are activated in response to many 
different abiotic and biotic stresses [38,39] and play a puta- 
tive role in plant defence responses to pathogens [17,40,41]. 
Indeed, Mai d 1 proteins are also known as pathogenesis- 
related proteins of class 10 (PR-lOs) [41]. However, the 
specific Mai d 1 genes involved in stress responses have 
not been identified and whether these genes are the same 
genes involved in allergy remains unknown. Thus, it is 
critical to determine whether breeding for low allergenicity 
would have consequences in terms of plant susceptibility 
against pathogens or abiotic stresses and whether high 
resistance to stresses would have consequences in terms 
of allergenicity. In addition, other stresses associated 
with growth and storage conditions affect the Mai d 1 
content in apples [2,3,42]. Therefore, exploiting our set 
of primers to study the variation of gene expression at 
different pre-/post-harvest conditions or after exposure 
to biotic/abiotic stresses might contribute to defining a 
"hypoallergenic protocol" for apple production to favour 
the reduction of symptoms in patients with apple allergies. 

Conclusions 

To date, patients with apple allergies deal with these 
condition through avoidance. To facilitate the normal 
consumption of essential components of a healthy diet 
containing apples, it will be useful to produce 'allergy- 
friendly' fruit. To obtain this goal, knowledge of the 
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identification of Mai d 1 isoallergens that cause allergy 
is needed. The qRT-PCR assay described in this work 
facilitates the examination of individual genes. The 
first application of this qRT-PCR assay showed that 
only a portion of the Mai d 1 genes is expressed in 
fruit. Moreover, the expressed genes showed great 
variation in expression among different tissues and 
cultivars. The data indicate that the presence and 
amount of specific Mai d 1 isoallergens determines 
allergenicity rather than total Mai d 1 content. For genes of 
specific interest, this assay might be further developed for 
allele (variant)-specific primer pairs, which, in turn, 
will promote the breeding of hypo-allergenic apple 
genotypes and support specific recommendations for 
fruit consumption, thereby reducing the impact of fruit 
allergies in a patients life. 

Methods 

Retrieval of Mai d 1 sequences and their alignment 

For the development of gene-specific primer pairs, all 
available sequence information for Mai d 1 was used at 
both the gene and allele levels. Coding DNA sequences 
(cds) were retrieved from the literature [6-8,17,19], 
the 'Golden Delicious' genome sequence [22,43] and a 
BLASTN search in the NCBI database [44] using keywords 
and Mai d 1 sequences as inputs. The coding sequences of 
31 Mai d 1 genes from the genome sequence of 'Golden 
Delicious' [6] were aligned using the Lasergene v8.0, 
MegAlign package (DNASTAR, Madison, WI, USA) and 
manually adjusted where necessary. A phenetic tree was 
produced using the same software, with default parameters 
and a neighbour-joining cluster algorithm (NJ). Single 
Nucleotide Polymorphisms (SNPs) among the different 
Mai d 1 sequences were identified. Subsequently, a 
complete alignment of all available Mai d 1 sequences 
was generated, and gene- and allele-differentiating SNPs 
were identified. 

Gene-specific primer design and validation 

Gene-specific primer pairs were designed using the 
software Primer3 [45]; the targeted gene-differentiating 
SNPs were located at the 3' end of at least one of the 
two primers (see Table 1) to ensure the primers ability 
to specifically amplify the target gene and, particularly, 
all known alleles of that target locus. Other criteria in 
the primer design included a primer length of 18-24 
nucleotides, a guanine-cytosine content of 20-80% and a 
RT-PCR amplicon length of 80-200 base pairs. Each 
primer pair was tested for the formation of homo- and 
heterodimers using the software PrimerSelect of Lasergene 
v8.0 (default settings). 

The gene specificity of each primer pair was validated 
in four ways. Firstly, for the in silico analysis, the primer 
sequences were blasted against the reference 'Golden 



Delicious' apple genome and the NCBI database. 
Secondly, to validate the primer specificity, the presence 
of a single PCR product after end-point PCR on genomic 
DNA was assessed. The PCR reactions were performed in 
a 17.5 \A volume containing 10 ng of DNA, 100 nM of 
gene-specific primers, IX reaction buffer, 1.5 mM MgCl 2 , 
100 [iM dNTPs and 0.5 Unit AmpliTaq Gold® DNA 
Polymerase (Applied Biosystems, Foster City, CA, USA). 
The reaction included an initial 10 min denaturation 
step at 95°C, followed by 33 PCR cycles (45 sec at 60°C, 
2 min at 72°C, and 30 sec at 95°C), with a final extension 
of 7 min at 72°C. The amplicons were visualised on an 
Image Station 440 CF (Kodak, Rochester, N.Y., USA) 
after electrophoresis on 1.5% (w/v) agarose gels and 
ethidium bromide staining. When non-specific bands 
were detected, the PCR conditions (primers concentration 
and annealing temperature) were optimised to obtain a 
single band. Thirdly, the primer specificity was validated 
through direct sequencing of all the amplicons from the 
gDNA of a set of 10 cultivars: 'Fiorina! 'Gala! 'Santanal 
'Elstar', 'Elise', 'Golden Delicious', 'Prima', 'Jonathan', 'Cox' 
and 'Ingrid Marie'. The sequencing reactions were 
performed at Bio-Fab Research, Pomezia, Italy and 
Greenomics, Wageningen, The Netherlands. The samples 
included tissues from low ('Elise' and 'Santana') and high 
('Golden Delicious') allergenic apple cultivars [2,3,46], that 
have been used in other studies of apple allergenicity and 
in studies on the effects of genetic modification on the 
expression of allergen genes ('Fiorina' and 'Gala') [24,25]. 
Cultivars currently used in apple breeding ('Cox', 'Elstar; 
'Golden Delicious', 'Jonathan' and 'Prima') were also 
included in these assessments. Subsequently, these 
sequences were analysed using Chromas Lite v.2.01, 
BLAST and MegAlign. Fourthly, the primer specificity 
was validated through qRT-PCR, using the shape of 
the melting curve as a criterion. Due to variations in 
length and nucleotide composition, each unique product 
was expected to have a unique melting temperature (Tm) 
and, consequently, a unique melting (or dissociation) curve. 
If the primer pair produced a single amplicon, the plot of 
the first derivative of the melting curve would contain 
a single sharp peak. 

Plant material, RNA extraction and cDNA synthesis 

The developed qRT-PCR assay was applied to the fruits 
of two apple cultivars ('Fiorina' and 'Gala'). The fruits 
were collected at the Cadriano Experimental Orchard of 
the University of Bologna (Italy) at the physiological 
ripening stage. For each genotype peel and flesh from 5 
apples were separately pooled, frozen in liquid nitrogen 
and stored at -80°C until RNA extraction. In particular, 
the peel was carefully separated from the flesh, leaving only 
a thin layer of flesh attached to the peel. Two different 
RNA extractions (biological replicates) were performed 
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from each pool. The cDNA was synthesised according 
to previously published methods [19]. 

Absolute quantification of the expression levels of the 
Mai d 1 genes 

Expression of Mai d 1 genes was characterised through 
qRT-PCR using a StepOne Plus Real-Time PCR instru- 
ment (Applied Biosystems) with a SYBR green-based 
assay. Similar to all other intercalating dyes, SYBR green 
binds to any double-stranded DNA; therefore, the speci- 
ficity of the primers was carefully assessed. Each reaction 
was performed in a total volume of 10 ul, containing 5 
ul of Power SYBR® Green Master Mix 2X, 70-100 nM 
of each primer, 3 ul of a 1:9 dilution of the cDNA and 
PCR-grade water. The reactions were incubated at 50°C 
for 2 min and at 95°C for 5 min, followed by 40 cycles of 
95°C for 15 sec and 60/63°C for 1 min, with data collection 
at each annealing step. The reactions were performed in 
triplicate (technical replicates). To ensure the specificity of 
the amplifications, each amplification reaction was followed 
by a melting phase, according to the default settings of the 
Step One Plus instrument (from 60°C to 95°C) and each 
melting curve was assessed for the presence of a single 
peak. The qRT-PCR assay included a standard curve for 
each target gene in the plate. Each standard curve com- 
prised 6 serial 1:10 dilutions (in duplicate) of the amplicons 
obtained from a fixed amount of gDNA using gene-specific 
primers, starting from 100 ng. The curve correlated fluores- 
cence signal (expressed as Ct values) to known amounts of 
amplicons in the 6 serial dilution samples and a regression 
line was plotted. Each assay also included a negative control 
performed in duplicate. The shape of the melting curves 
obtained from the standard curve dilutions (gDNA) and 
the samples (cDNA) indicates whether the amplified 
products are homogeneous and the melting temperature 
provides confirmation that the correct product has been 
specifically amplified. 

The amplification efficiency (E) of each primer pair 
was estimated using the slope of the regression line, 
according to the equation: E=10 ( " 1/slope) - 1 [47]. The 
relative PCR efficiencies were evaluated for each target 
gene in relation to the reference gene to select an appro- 
priate method for the analysis of the raw qRT-PCR data. 
To obtain the gene expression results, the raw data was 
transformed using the "Standard curve method" [23] 
and reported as relative expression levels; the transcript 
levels of the Mai d 1 genes were normalised to the 
transcript levels of actin (primer sequences from [48]) and 
expressed as Arbitrary Units (A.U.). Actin was chosen after 
the comparison of its gene expression to other two putative 
reference genes: (1) UBC (UBiquitin-Conjugating enzyme, 
MDP0000223660), forward primer 5'-CGAATTTGTCCG 
AAGGCGT-3'; reverse primer: 5'-CAATGATTGTCACAG 
CAGCCA-3', and (2) GAPDH (Glyceraldehyde 3-phosphate 



dehydrogenase, MDP0000835914), forward primer 5'-ATT 
GGCAGTGTGCGACGTT-3' and reverse primer 5'-GGA 
GGAGTCAATGGTGGAGGA-30. The comparison of the 
three putative reference genes was done in two biological 
replicates of peel and flesh of 'Fiorina and 'Gala . The mean 
normalised expression levels and the standard error of the 
mean (SEM) were calculated using the two biological 
replicates. Univariate analyses of the differences between 
the mean values were performed using analysis of variance 
(ANOVA) with a 0.05 significance level. 

Additional files 



Additional file 1: List of Mai d 1 sequences included in the 
alignment. 

Additional file 2: Alignment of the 31 coding sequences of the Mai 
d 1 isoallergen genes. Alignment of the 31 coding sequences of the 
Mai d 1 isoallergen genes retrieved from the 'Golden Delicious' genome 
sequence. The alignment was performed using MegAlign (DNASTAR 
Lasergene v8.0). Each sequence was reported using the name of the 
related gene and the accession number from the Apple GBrowse - Malus 
x domestica v1.0 [22], The mismatched residues in the consensus 
sequence are highlighted in white. 

Additional file 3: Results of direct sequencing. The results of the 
BLASTN analysis of the sequences retrieved from the direct sequencing 
of Mai d 1 amplicons from 10 different apple cultivars. The table lists the 
targeted gene, the name of the primer pair and the ID of the most 
similar sequence in the database or in the Gbrowse of the apple genome 
for each of the cultivars [22]; the presence/absence of the SNPs in these 
reference sequences is also indicated. 
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